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SECTION  1 


INTRODUCTION 

The  active  and  downstream  plasma  processing  of  semiconductor  devices  and  surface 
modification  of  aerospace  materials  is  fundamental  to  Air  Force  and  other  DOD  applications. 
There  is  a  great  need  for  a  more  fundamental  understanding  of  these  sorts  of  cold,  reactive 
plasmas.  There  is  also  a  need  for  the  development  of  practical  diagnostics  for  plasmas  actually 
used  or  potentially  of  use  in  materials  processing.  Of  these  techniques,  laser  diagnostics  are 
almost  certainly  among  the  most  powerful  and  promising. 

A  particularly  striking  confluence  of  fundamental  and  applied  interests  occurs  with  the 
use  of  H2  plasmas  in  the  processing  of  UI-V  electronic  devices.  On  the  fundamental  side,  there 
is  probably  no  system  with  simpler  atomic  and  molecular  species.  On  the  other  hand,  the  H; 
system  exhibits  all  the  properties  of  typical  reactive  plasmas.  The  obtainment  of  a  good 
fimdamental  understanding  of  the  Hj  RF  plasma  would  seem  a  necessary  prelude  to  the  complete 
understanding  of  more  complex  plasma  systems. 

However,  just  as  importantly,  the  H2  plasma  has  immediate  application  in  the  processing 
of  111-V  materials.  Thus,  there  is  every  reason  to  believe  that  these  same  diagnostic  techniques 
can  be  applied  to  elucidate  the  mechanisms  of  the  processing  of  III-V  materials  in  H2  plasmas 
or  other  plasmas  containing  hydrogenic  species.  Ultimately,  this  diagnostic  may  provide  a 
unique  approach  to  understanding  the  plasma  induced  surface  chemistry  in  these  systems.  While 
a  basic  knowledge  of  the  plasma  processes  is  a  necessary  prerequisite  for  a  complete 
understanding  of  these  systems,  it  should  be  possible  to  make  concurrent  progress  in  controlling 
aixl  manipulating  plasma  processing  by  applying  laser  diagnostics  while  these  diagnostics 
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themselves  are  in  development. 

Hj  and  H2/CH4  plasmas  are  used  to  directly  etch  GaAs  and  also  InP.'  Of  at  least 
comparable  interest  is  the  use  of  H2  plasmas  to  passivate  devices.  For  example,  while  other 
etchant  gases,  e.g.  CI2,  appear  at  the  momem  preferable  for  etching  some  IIl-V  materials,  e  g. 
GaAlAs,  passivation  in  an  H2  plasma  of  surface  defects  created  in  other  etching  processes,  seems 
very  effective.^  Moreover,  device  structures  can  be  produced  by  selective  (using  masks)  H2 
passivation  of  dopants.  Examples  include  the  production  of  a  variety  of  GaAs  FET  with  a 
hydrogenated  channel  and  stripe  geometry  AlGaAs-GaAs  quantum  well  lasers  using  masked 
hydrogenation  to  produce  current  guiding  through  the  heterostructure.* 

Unintentional,  sometimes  deleterious,  hydrogenation  can  also  occur  in  processing 
environments.  For  example,  H2  fluorocarbon  mixtures  are  sometimes  used  for  etching  SiOj  and 
SiN,  layers  on  III-V  semiconductors.^  Even  "accidents”  such  as  introduction  of  moisture  into 
a  plasma  can  lead  to  H-atom  production. 

The  report  describes  the  systematic  development  of  a  reliable  and  accurate  diagnostic  for 
probing  of  the  most  abundant  reactive  species  in  a  cold,  processing  hydrogen  plasma,  the  H- 
atom.  The  diagnostic  provides  a  spatially  ( » 100/x)  and  temporally  ( » 10  nsec)  resolved  relative 
( a  10  percent  precision)  and  absolute  ( « 2S  percent  precision)  measure  of  the  H-atom 
concentration  within  the  plasma.  This  diagnostic  is  based  upon  the  two-photon  absorption  laser- 
induced  fluorescence  (TAUF)  teclmique.  It  is  described  in  detail  in  the  second  section  of  this 
report. 

As  we  have  developed  the  diagnostic,  we  have  applied  it  to  an  elucidation  of  the  cold  H2 
plasmas  typical  of  processing  environments.  In  Section  3  A,  we  describe  the  general 
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dependences  of  the  H-atom  concentration  on  discharge  parameters.  In  Section  3.B,  we  report 
H-atom  spatial  profiles  in  continuous  and  pulsed  H2  discharges.  We  further  present  a  model  of 
the  H-atom  production  and  concentration  evolution  that  is  in  strikingly  good  agreement  with  the 
measured  profiles.  We  then  further  simulate  a  processing  environment  by  introducing  GaAs  and 
Si  wafers  into  the  reactor  and  observing  the  modifications  caused  thereby  to  the  H-atom 
concentrations.  The  H-atom  concentration  model  for  an  unloaded  reactor  is  extended  and 
successfully  describes  the  situation  where  wafers  are  present  on  the  electrode.  In  Section  3.D, 
we  describe  how  the  TALIF  linewidths  can  be  used  as  a  sensitive  probe  of  temperature  and  gas 
purity  within  the  reactor. 

Section  4  is  a  short  surtunary  of  conclusions  reached  and  suggestions  for  further  work. 
The  Appendix  lists  publications  and  presentations  based  upon  work  supported  by  this  contract. 
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EXmRIMENTAL  DETAILS 
A.  LASER  CONnCURATION 

A  schematic  of  the  general  experimental  scheme  used  for  the  TALIF  H-atom  detection 
is  shown  in  Fig.  1.  The  schematic  layout  shows  an  RF  reactor  cell  power  supply  system,  the 
laser  system,  and  the  associated  acquisition  and  control  electronics.  The  laser  light  generation 
configuration  will  now  be  described  and  the  acquisition  electronics  are  discussed  below. 

The  TALIF  transition  for  atomic  hydrogen  is  from  the  Is^S,/!  ground  state  to  the 
3d^D5,2  3,2  and  3s^S„2  excited  states  at  97492  cm Beth  excited  state  levels  are  populated  due  to 
the  energy  level  degeneracy  in  hydrogen  for  different  angular  momentum  states  and  due  to  the 
two-photon  absorption  selection  rule  of  AL  =  0,  ±2.*'  The  transition  requires  two  205.14-nm 
photons  to  excite  the  atomic  ground  state.  The  TALIF  transition  has  also  been  observed  in 
deuterium  at  97519  cm  ‘  in  the  course  of  this  research. 

The  excited  state  is  only  weakly  populated  by  the  discharge  impact  kinetics  due  to  a  rapid 
(  =  15  nsec)  radiative  lifetime.  Thus,  the  TALIF  excitation  can  be  expected  to  effectively  alter 
the  population  of  the  upper  state  and  its  emission,  although  significant  photon  density  is  required 
given  the  small  two-photon  absorption  cross  section. 

TALIF  measurements  for  H-atoms  can  be  performed  using  a  high  pressure  H2-filIed 
Raman  cell  for  frequency  shifting  or  by  sum-frequency  mixing  (SFM)  in  a  beta-barium  borate 
(BBO)  crystal.  For  our  purposes  the  use  of  SFM  with  a  BBO  crystal  is  preferred  since  SFM 
generates  far  more  than  205-nm  light.  The  light  generation  scheme  depicted  in  Fig.  1  was  used 
in  all  but  our  initial  experiments.  The  probe  laser  is  generated  by  pumping  a  Quanta-Ray  PDL-2 
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Figure  1.  Schematic  for  production  of  20S-nm  light  for  H-atom  TALIF  excitation  and  * 

detection.  Note  that  the  photomultiplier  (PMT)  is  looking  perpendicular  to  the 
205-nm  beam  even  though  it  is  depicted  in  the  plane  of  the  Figure. 
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dye  laser  with  the  second-harmonic  532-nm  output  of  a  Quanta-Ray  DCR-2A  Nd.  YAG  laser 
operating  at  20  Hz.  The  dye  laser  operates  with  Exciton  Sulfoiiiodamine  640  dye  at  615  nm.  The 
dye  laser  output  is  frequency  doubled  using  a  Quanta-Ray  WEX  KDP  crystal  providing  an  output 
of  7  mJ/pulse  at  307  nm  and  a  residual  14  mJ/pulse  at  615  nm.  Upon  exiting  the  WEX  system, 
the  polarization  of  the  307-nm  beam  is  rotated  to  coincide  with  the  polarization  of  the  615-nm 
beam.  Both  beams  are  collimated  with  a  3:1  telescope  and  combined  via  a  dichroic  mirror.  The 
collinear,  collimated  beams  are  input  to  the  BBO  crystal  manually  angle  tuned  to  provide  SFG 
at  205  nm.  The  spot  size  diameter  of  the  input  beams  is  3  mm.  The  SFG  output  from  the  BBO 
crystal  can  be  as  high  as  800  microJoules/pulse  with  careful  alignment  and  overlap  of  the  two 
input  beams.  The  output  pulse  energy  is  carefully  measured  with  a  calibrated  Molectron  Model 
J3-05  Joulemeter.  Using  the  BBO  SFM  technique  gives  over  10  times  the  output  of  the  Hj 
Raman  shift  technique.  The  conversion  efficiency  of  BBO  at  this  wavelength  has  been  quoted 
possibly  as  high  as  12  percent  in  the  literature’,  and  the  efficiency  of  our  system  was  measured 
at  1 1  percent  which  is  predominantly  dictated  by  the  input  beams  profile  quality  and  alignment. 
Upon  leaving  the  BBO  crystal  the  residual  input  pump  beams  are  dispersed  by  a  Pellin-Broca 
prism  and  beam  dumped.  The  TALIF  probe  beam  is  tightly  focused  by  a  quartz  lens  through 
a  Supers!  1  1  quartz  window  into  the  discharge  reactor  of  interest. 

B.  SIGNAL  DETECTION  AND  ANALYSIS 

Refinement  and  application  of  the  H-atom  TALIF  diagnostic  was  performed  using  a 
reactor  that  has  a  variable  gap  cylindrical  planar  electrode  design  with  water-cooled  stainless 
steel  or  aluminum  electrodes.  The  discharge  load  power  can  be  measured  with  a  high  voltage 
probe  and  a  current  loop  detector  whose  output  is  acquired  on  a  digital  oscilloscope  and 
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multiplication  averaged  upon  waveform  transfer  to  a  personal  computer.  The  reactor  cell  is 
mounted  on  a  stepper  motor  positioned  X-Z  translation  stage.  See  Fig.  1,  X  is  the  in-plane 
direction  perpendicular  to  the  laser  beam  and  Z  is  perpendicular  to  the  plane  of  the  figure. 
This  configuration  allows  spatial  profile  measurements  as  the  cell  is  translated  relative  to  the 
fixed  laser  probe  focal  point  as  viewed  by  the  fluorescence  collection  optics. 

The  excited  state  fluoresceiKe  from  the  3d^D,^  j;,  and  levels  to  the  on 

the  Balmer  alpha  line  (6S6.3  nm)  transition  is  detected  perpendicular  to  the  primary  discharge 
cell  axis  and  the  laser  probe  beam  with  focused  collection  optics.  The  TALIF  passes  through 
a  lO-nm  FWHM  interference  filter  centered  at  656  nm  and  focused  on  the  photocathode  of  an 
ITT  Model  4123  gated  photomultiplier  tube  (GPMT).  Detection  of  H-atom  TALIF  is  also 
easily  achieved  (provided  the  average  DC  current  does  not  exceed  the  rating  of  the  particular 
photomultiplier  tube)  by  using  continuous  (CW)  photomultipliers  such  as  the  RCA  1P28  or 
Thom-EMI  9659  tubes. 

The  signal  collected  by  the  either  a  GPMT  or  CW  photomultiplier  is  fed  directly  to  the 
input  Tektronix  DSA  601  Digital  Signal  Analyzer  or  other  similar  high-speed  digital 
oscilloscope.  The  data  acquisition,  dye  laser  tuning,  and  experiment  timing  is  controlled  by  a 
Stanford  research  System  Digital  Delay  Generator  in  combination  with  an  IBM-compatible  386 
personal  computer.  The  experiment  timing  sequence  allows  for  temporally-dependent 
measurements  of  H-atom  TALIF  over  all  relative  phases  of  the  RF  discharge  driving  cycle.  The 
use  of  a  GPMT  and  a  narrow-band  interference  filter  in  the  collection  optics  allows  detection 
of  the  TALIF  response  even  in  the  presence  of  the  bright  plasma-induced  emission  (PIE).  In 
working  with  the  larger  reactor  systems  it  was  found  that  even  using  a  CW  photomultiplier  was 
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sufficient  provided  the  anode  current  of  the  tube  was  carefully  monitored  and  did  not  exceed 
the  maximum  rating. 

A  representative  H-atom  TALIF  signal  trace  is  shown  in  Figure  2.  The  TALIF 
resonance  transition  is  shown  for  a  typical  set  of  discharge  operating  conditions.  As  shown,  the 
TALIF  signal  has  an  excellent  signal/noise  ratio  wdiich  is  typical  of  the  TALIF  signal  over  all 
phases  of  the  RF  cycle  relative  to  the  PIE.  The  full-width-half-maximum  (FWHM)  two-photon 
linewidth  of  the  TALIF  resonance  is  typically  3.5  ±  1  cm  '.  The  dye  laser  linewidth  at  615  nm 
is  approximately  0.6  cm as  measured  independently  by  optogalvanic  detection  of  neon  and 
iron  transitions  in  a  hollow-cathode  glow  discharge  lamp.  Multiplying  this  by  a  factor  of  6 
(third  harmonic  of  615  nm  times  two  photons)  accounts  for  the  FWHM  linewidth  of  the  TALIF 
transition. 

Note  the  associated  Doppler  linewidth  is  <1.0  cm  '  for  the  two-photon  absorption  in 
hydrogen  at  a  neutral  atom  temperature  of  300-500  Kelvin  which  might  be  expected  in  the 
center  of  the  discharge  gap  at  zero  degrees  phase  angle  in  the  RF  current  cycle.  Thus  our 
TALIF  measurements  do  not  have  a  sufficiently  small  laser  probe  linewidth  to  measure  H-atom 
linewidths  with  typical  thermal  energies  without  etalon  narrowing  of  the  dye  laser  fundamental 
output. 

In  some  experiments  (see  Section  3.D)  the  linewidth  of  the  dye  laser  was  narrowed  by 
an  intracavity  etalon  to  <«0.06  cm  '.  The  dye  laser  is  scanned  by  pressure  tuning,  wliich  entails 
evacuating  and  pressurizing  with  dry  Nj  a  box  surrounding  the  etalon  and  grating.  This  system 
is  controlled  by  a  solenoid  valve  switched  by  a  relay  that  is  actuated  by  the  data  acquisition 
software. 
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c. 


H-ATOM  CONCENTRATION  CALIBRATION  PROCEDURE 


Having  developed  a  laser  diagnostic  to  detect  H-atoms,  it  would  be  useful  to  place  the 
amount  of  atoms  detected  on  an  absolutely  calibrated  scale  and  determine  actual  local 
concentrations  of  atoms  directly  in  the  plasma.  By  directing  an  external  source  of  H-atoms  into 
a  central  position  in  the  volume  of  a  plasma  reactor,  the  TALIF  signal  derived  from  this  source 
can  be  made  to  simulate  detection  of  plasma-generated  atoms.  The  external  source  of  atoms, 
usually  generated  via  microwave  discharge,  provides  a  steady,  reproducible,  and  controllable 
TALIF  response.  Simultaneously,  small  quantities  of  molecules  known  to  efficiently  react  and 
combine  with  the  atomic  species  can  be  added,  and  a  corresponding  decrease  in  the  TALIF 
signal  intensity  can  be  observed  as  the  atoms  are  depleted.  This  reaction  or  titration  technique 
can  give  a  reliable  means  of  determining  the  concentration  of  atomic  species  corresponding  to 
a  given  amount  of  TALIF  signal. 

Calibration  of  the  fluorescence  signal  collected  by  the  PMT-filter  combination  was 
performed  using  in-situ  titration  of  the  H-atoms  by  chemical  reaction  with  NOj.  The  titration 
reaction  proceeds  as* 

+  NOj  -  Off  +  NO 

with  a  reaction  rate  constant  of  k29,  =  1.3  x  10  '“  cm’  s  '  w^cre  the  rate  constant  is  given  for 
room  temperature  (T=298K).  There  is  also  a  secondary  reaction  of  OH  with  Hj  given  by 


with  a  reaction  rate  constant  kj,,  =  1.6  x  lO  ”  cm^  s'.  However,  provided  the  ambient  Hj 
concentration  remains  below  1  Torr  (3.3  x  10'^  cm'^)  this  reaction  will  not  appreciably  perturb 
the  main  reaction.  Sufficient  mixing  time  must  be  allowed  to  ensure  the  reaction  goes  to 
completion.  To  have  a  reaction  go  to  99  percent  completion  implies  that  the  mixing  time,  t„,„„ 
must  be  greater  than  S  characteristic  reaction  times,  t.  In  other  words,  the  percentage  of 
complete  reaction  is 

l-e-<'’  =  l-e-*"'"  =  l-e-’=0.9932-99% 

The  characteristic  reaction  time  is  the  inverse  of  the  reaction  rate  constant  times  the 
concentration  of  reactants.  For  example,  in  the  case  of  the  reaction  in  Equation  1,  at  the  point 
of  complete  reaction  (here  called  endpoint)  a  measured  concentration  of  NO^  is  used  to  react 
completely  with  an  equal  concentration  of  H-atoms.  If  the  concentration  of  [NOj]  used  is  I  x 
10’’  cm then  the  required  mixing  time  must  be  at  least 

r^i5x=5(ik2^[A102])'‘=7.7  microseconds  U) 

In  the  titration  assembly  described  below,  the  main  gas  flow  rates  are  adjusted  to  reduce  the 
slug  flow  gas  velocity  in  the  mixing  zone  and  raise  the  mixing  time  to  approximately  4 
milliseconds,  clearly  satisfying  the  condition  in  Equation  2. 

The  titration  assembly  is  shown  in  Fig.  3.  At  the  side  inlet  a  gas  mixture  of  Hj/Ar  is 
passed  through  a  microwave  discharge  produced  by  an  air-cooled  Evenson  microwave 
discharge  cavity  which  dissociates  the  molecular  gas  mixture.  The  presence  of  argon  helps  to 
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titration  of  H-atoms. 


induce  more  complete  dissociation.  The  discharge  section  is  composed  of  quartz  for  both  heat 
resistance  and  to  minimize  reactions  of  the  microwave  plasma  with  the  walls.  Further 
downstream,  NO2  is  introduced  through  a  glass  coaxial  tube  which  passes  the  gas  to  a 
perforated  mixing  zone  approximately  2  cm  in  length  and  2  cm  in  diameter.  The  inner  tube  is 
covered  by  a  Teflon  sleeve  to  reduce  the  amount  of  recombination  of  H-atoms  on  the  tube 
surface.  The  assembly  is  placed  through  a  3/4  inch  diameter  (l.D.)  Cajon  0-Ring  fitting  welded 
to  a  2-3/4  inch  Conflat  copper  gasket  sealed  flange.  The  Conflat  flange  is  mated  to  a  port  on 
a  given  reactor  and  the  length  of  the  titration  assembly  allows  the  exit  from  the  mixing  zone 
to  be  placed  at  the  center  of  the  reactor  and  near  the  focal  point  of  the  probe  laser.  For  a 
titration  calibration  measuement,  a  steady-state  concentration  of  H-atom  is  produced  by  the 
microwave  cavity  discharge  source  and  NO^  is  introduced  at  known  mass  flow  rates.  The 
TALIF  fluorescence  is  monitored  and  recorded  manually  as  a  function  of  addition  of  NO^.  A 
typical  titration  data  curve  is  shown  in  Fig.  4. 

At  the  endpoint  (minimum  observable  TALIF  signal),  the  total  pressure  in  the  reactor 
is  recorded.  A  straight  line  is  least-squares  fit  to  the  first  several  titration  data  points,  and  the 
exact  NOj  endpoint  mass  flow  rate  is  extrapolated  numerically.  Combining  this  information, 
we  can  use  the  following  equation,  where  at  the  endpoint 


[H]^  =  tiVOjlo  =(3.3xlO“cm-’rorr-‘)*f^!^^^ 


•P 


(3) 


where  [H]o  is  the  concentration  of  H-atoms  titrated,  [NOj],,  is  the  concentration  of  NO^  used, 
FfNOj)  is  the  mass  flow  rate  of  NOj  measured  in  standard  cubic  centimeters  per  minute 
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Figure  4.  Absolute  calibration  curve  of  TAUF  intensity.  Inset:  typical  titration  data  show  i 

a  linear  decrease  of  the  TALIF  intensity  upon  addition  of  NO2  in  the  flow  reactor. 

The  NO2  concentration  at  the  jc-intetcept  equals  the  H-concentration  corresponding 
to  the  TAUF  intensity  at  the  y-intercept.  The  H-concentrations  are  then  plotted 
against  the  observed  UF  iidensity  in  the  main  plot  to  obtain  the  calibration  curve. 
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(seem),  IF  is  the  total  mass  flow  rates  of  all  the  gases  used  summed  together,  and  P  is  the  total 
pressure  measured  in  Torr.  The  numerical  coefficient  is  used  to  convert  Torr  into  number 
density,  per  cm^  Titration  curves  can  be  taken  for  different  initial  atomic  concentrations,  PMT 
gain  settings,  and  probe  laser  energy.  The  calibrations  curves  (see  Fig.  4)  correlate  an 
experimental  TALIF  signal  with  a  known  concentration.  Typically  the  calibrations  are 
reproducible  to  IS  percent  and  believed  absolutely  accurate  to  25  percent.  The  minimum 
detectable  concentration  is  extrapolated  to  be  10'^  cm'\ 
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H-ATOM  DIAGNOSTIC  RESULTS 

A.  h-atom  dependence  upon  discharge  parameters 

It  is  useful  to  get  a  general  overview  of  the  dependence  of  the  TALIF  signal  upon  laser 
and  discharge  parameters  to  characterize  its  usefulness  as  a  plasma  diagnostic.  Fig.  S  shows 
the  TALIF  signal  response  as  a  function  of  incident  205-nm  laser  probe  pulse  energy.  The 
TALIF  absorption  should  depend  quadratically  upon  the  laser  pulse  energy  for  a  fixed 
interaction  volume.  The  curve  through  the  data  points  in  Fig.  S  is  a  quadratic  polynomial 
regression  fit  to  the  data  and  shows  agreement  with  the  expected  TALIF  signal  response.  The 
lack  of  a  saturation  behavior  in  the  data  indicates  that  the  TALIF  signal  response  is  not 
affected  by  competing  processes  such  as  multiphoton  ionization.  In  addition,  linewidth 
measurements  at  zero  degrees  phase  show  no  power  broadening  effects  for  this  range  of  laser 
pulse  energies. 

The  H-atom  TALIF  as  a  function  of  RF  discharge  load  power  is  shown  in  Fig.  6.  The 
TALIF  signal  has  been  acquired  for  both  100%  Hj  and  Dj  gases  at  1  Torr,  50  seem,  100  kHz 
RF,  and  zero  degrees  RF  current  phase.  The  signals  shown  in  Fig.  6  are  corrected  for 
differences  in  laser  pulse  energies  using  the  data  shown  in  Fig.  5.  Interestingly,  the  deuterium 
TALIF  at  97519  cm  '  shows  a  factor  of  2  decrease  relative  to  the  hydrogen  TALIF  owing 
either  to  differences  in  the  Dj  versus  dissociation  and  recombination  rates  or  the  two-photon 
absorption  cross  sections.  However,  both  TALIF  signals  appear  to  saturate  at  about  75-W  RF 
power,  perhaps  due  to  a  saturation  of  electron-impact  dissociation  formation  of  ground  state 
atoms. 
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TALIF  SIGNAL  RF  POWER  DEPENDENCE 
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The  dependence  of  the  TALIF  on  discharge  pressure  for  hydrogen  and  deuterium  is 
demonstrated  in  Fig.  7.  Here  the  RF  discharge  load  power  was  kept  fixed  and  both  TALIF 
signals  normalized  to  the  same  20S  nm  laser  pulse  energy.  The  TALIF  has  a  rapid  increase 
from  0.2  Torr  to  a  maximum  at  about  1.75  Torr  followed  by  a  slow  decline  through  10  Torr.  ^ 

Both  the  H2  and  D2  TALIF  show  comparable  behavior  with  the  signal  response  following  the 
maximum  being  interpreted  as  a  loss  of  fluorescence  due  to  excited  state  quenching  and/or  an 
increase  in  atom  recombination  losses  as  the  collision  frequency  increases  with  the  higher  * 

pressures. 

The  phase  dependence  of  the  H  atom  concentration  at  100  kHz  is  shown  in  Fig.  8.  The 
data  were  taken  in  the  middle  of  the  electrode  gap  at  1.0  Torr  of  50  percent  hydrogen-50 
percent  helium  and  6-W  load  power.  The  upper  trace  is  the  TALIF  signal;  the  response  is 
constant  within  experimental  error  over  all  phase  angles.  No  phase  dependence  was  observed  • 

even  for  frequencies  as  low  as  20  kHz;  this  implies  the  H  atom  ground-state  lifetime  in  these 
discharges  is  greater  than  50  psec.  The  lower  trace  is  the  relative  response  of  the  PIE  over  the 

» 

cycle.  A  crucial  point  in  making  this  measurement  is  the  necessity  of  ensuring  that  the  GPMT 

does  not  saturate  at  phases  when  the  PIE  is  very  bright.  This  means  that  the  GPMT  must  be 

run  at  a  low  gain  (10^-10^)  in  order  to  achieve  a  linear  response  from  the  GPMT.  Indeed,  in  B 

Fig.  8,  at  just  about  the  level  of  experimental  uncertainty,  there  appears  to  be  a  very  small  dip 

in  the  TALIF  signal  at  the  maximum  of  the  PIE,  presumably  due  to  an  incipient  saturation 

» 

effect. 

The  H-atom  TALIF  signal  versus  H2  concentration  in  helium  is  shown  in  Fig.  9.  The 
TALIF  signal  is  nonlinear  over  most  of  the  Hj  percentage  concentrations.  However,  the  inset  ^ 
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Figure  9.  H-atom  TALIF  signal  versus  H,  coaceniration  in  helium.  The  total  discharge  cell  pressure  is  1  Torr.  The  inset  shows 
an  expanded  portion  of  the  plot  for  concentrations  below  10%. 


shows  that  below  «1  percent  the  TALIF  response  is  nearly  linear  and  passes  through  the  origin. 
B.  H-ATOM  PROFILES  IN  UNLOADED  DISCHARGES 

i)  Observations  in  continuous  and  pubed  discharges 

A  typical  axial  profile  of  H-atom  concentration  at  different  positions  across  the 
interelectrode  space  of  a  continuous  RF  discharge  and  above  the  center  of  the  electrode  is 
shown  in  Fig.  10.  It  was  taken  with  the  laser  frequency  sitting  at  the  peak  of  the  transition 
while  the  translation  stage  is  stepped  to  different  positions  across  the  I9-min  gap.  All  the 
spatial  profiles  shown  in  this  work  are  corrected  for  the  occultation  of  the  detection  solid  angle 
from  the  electrode  surface  as  the  distance  between  the  laser  beam  and  the  electrode  surface  gets 
smaller  than  3.S  mm.  A  calculation  of  the  reduction  of  the  solid  angle  from  which  we  collect 
the  6S6-nm  fluorescence  is  checked  by  experiment.  We  fill  the  reactor  cell  with  CjHj  at  a 
pressure  of  SO  mTorr  and  we  detect  H-atoms  produced  by  photolysis  of  the  acetylene  molecule 
with  205-nm  photons.’  Since  the  detection  of  H-atoms  takes  place  immediately  after  their 
production,  within  a  laser  pulse  duration  (10  ns),  the  concentration  is  independent  of  the 
position  of  production  anywhere  in  the  interelectrode  space.  The  reduction  in  the  TALIF  signal 
that  we  observed  is  a  result  of  the  viewing  reduction  of  the  solid  angle  as  the  laser  beam 
approaches  the  electrode  surface,  which  is  determined  by  the  detection/collection  geometry. 
We  have  found  that  the  experimental  determination  of  the  geometrical  view  factor  agrees 
within  10  percent  with  that  calculated  by  geometrical  considerations. As  shown  in  Fig.  10, 
a  concentration  of  about  8x10'^  atoms/cm^  stays  relatively  constant  over  a  wide  region  around 
the  center  of  the  discharge.  As  we  approach  the  electrode  surfaces,  however,  the  concentration 
drops  sharply  to  zero  within  3-4  mm  away  from  each  electrode  surface.  We  discuss  below  the 
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mechanism  that  accounts  for  the  observed  concentration  characteristics,  where  we  model  the 
temporal  evolution  of  the  concentration  i^files.  The  solid  line  in  Fig.  10  results  from  the 
proposed  model  (see  3.B.2)  and  shows  a  reasonably  good  agreement  with  the  measurements. 

Spatial  profiles  of  hydrogen  atoms  are  also  taken  in  a  pulsed  discharge  at  different  times 
after  the  discharge  is  turned  on.  In  Fig.  1 1,  we  show  the  spatial  distribution  of  Hatoms  at  0.2S 
ms,  1  ms,  and  2.S  ms.  The  figure  describes  the  temporal  evolution  of  H-atoms.  The  profile 
at  0.2S  ms  clearly  shows  that  H-atoms  are  generated  at  regions  centered  at  about  3.4  mm  away 
from  each  electrode  surface.  The  1.0-  and  2.5-ms  traces  show  how  the  steady-state  distribution 
builds  up.  H-atoms  exit  the  production  regions  and  fill  the  interelectrode  space  until  they 
approach  the  metallic  electrode  surfaces,  where  they  stick  or  recombine.  The  profile  at  2.5  ms 
closely  resembles  that  shown  in  Fig.  10,  taken  in  a  continuous  RF  discharge,  indicating  that 
a  steady  state  is  essentially  established  in  the  discharge  at  this  point  in  time. 

The  profiles  in  Fig.  12  represent  typical  H-atom  distributions  after  the  discharge  is 
turned  off  at  tg.  The  spatial  profiles  shown  are  taken  at  elapsed  time  t  after  tg,  with  t=0.25  ms, 
t=0.5  ms  and  t=l  ms.  It  is  clear  that  the  electrodes  represent  a  significant  sink  for  H-atoms  in 
the  system.  The  solid  lines  shown  in  Figs.  II  and  12  represent  best  fitted  theoretical 
concentration  profiles  which  we  discuss  below. 

Fig.  13  shows  a  temporal  profile  of  the  hydrogen  atom  concentration  at  the  center  of 
the  discharge.  The  discharge  is  turned  off  at  to=2  ms,  a  time  close  to  the  establishment  of  an 
equilibrium  concentration  in  the  discharge.  The  decay  of  hydrogen  atoms  follows  a  single 
exponential  law  as  indicated  by  the  linearity  of  the  logarithmic  fit  displayed  also  in  Fig.  13, 
that  is  independent  of  the  atomic  concentration  at  the  time  the  discharge  was  turned  off.  The 
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Figure  1 1 .  Temporal  evolution  of  H-atcmi  spatial  profiles  in  a  pulsed  discharge;  experimental 
data  points  and  simulated  results  (solid  lines)  at  t=0.2S,  1.0,  and  2.5  ms.  The 
pressure  was  nuintained  at  3  Torr  and  the  power  was  5  W. 

i 


26 


i 


Time  (ms) 

Figure  13.  Temporal  profile  of  H-alom  concentratioii  taken  in  the  tftei^low  of  a  2-in8  pulse 
discharge  at  the  center  of  the  interelectrode  iqwce.  A  single  exponential  (as 
shown  by  log  plot)  can  fit  well  the  decay  of  H-concentration  at  times  S200^ 
after  the  discharge  pulse  ends.  Discharge  conditions  as  in  Fig.  10. 


exponential  decay  reveals  that  first  order  kinetics  are  re^nsible  for  the  decay  of  hydrogen 
atom  concentration  after  the  dischai^e  is  turned  off,  widi  a  time  constant  t=1.03±0.08  ms. 
Therefore,  gas  phase  recombination  of  the  hydrogen  atoms  is  not  significant  since  it  would 
result  in  second  order  kinetics  with  respect  to  H-atom  concentration. 

Finally,  we  replace  deuterium  for  hydrogen  to  check  for  any  quantitative  differences  in 
behavior  due  to  their  mass  difference.  From  the  decay  of  the  atomic  deuterium  concentration, 
we  determine  a  time  constant  t=1.66±0.10  ms.  The  comparison  of  the  time  constants  of 
hydrogen  and  deuterium  indicates  an  £q)proximate  m'*'^  dependence  on  the  mass  of  the 
decaying  p>articles. 

ii)  Model 

The  experimental  observations  presented  in  the  previous  section  guided  us  in  the 
formulation  of  a  model  that  incorporates  the  following  assumptions;  First,  the  gas  phase 
recombination  is  negligible  as  indicated  by  the  simple  exponential  decay  of  the  atomic 
hydrogen  population.  Second,  the  H-atoms  are  destroyed  only  on  the  metallic  electrode 
surfaces  which  induce  catalytical  recombination  and  therefore  represent  sinks  for  H-atoms. 
Third,  neutral  diffusion  governs  the  movement  of  the  H-atoms  towards  regions  of  reduced 
atomic  concentration  in  both  the  discharge  and  the  afterglow.  Finally,  the  production  of  H- 
atoms  takes  place  in  narrow  regions  in  the  discharge.  These  regions  roughly  coincide  with  the 
visibly  brightest  portion  of  the  plasma,  indicative  of  regions  of  high  excitation  and  dissociation. 

A  quantitative  formulation  of  our  model  is  given  as  follows.  Let  x  be  the  coordinate 
on  the  axis  perpendicular  to  the  electrode  surfaces  (sec  Fig.  1).  Because  of  the  symmetry  of 
our  discharge,  we  assume  a  uniform  spatial  distribution  in  y-  and  z-  directions,  so  that  the 
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atomic  hydrogen  concentration  p  depends  only  on  x  and  t.  We  approximate  the  source  regions 
with  delta  functions  centered  at  a  distance  \  away  from  each  electrode  surface.  Clearly,  this 
assumption  is  oversimplified  but,  as  we  shall  see,  it  yields  reasonable  results.  Thus,  we  solve 
the  one-dimensional  diffusion  equation:" 


dp/dt  -  D  ^pldx^  =  Po  (  5(x-xJ  +8(x-a+Xo)) 


where  D  is  the  diffusion  coefficient  for  the  atomic  species,  Pg  is  the  production  rate  of  the 
sources  and  a  is  the  interelectrode  separation.  We  will  assume  for  simplicity  that  the  electrodes 
are  perfectly  absorbing  walls  for  the  hydrogen  atoms,  and  we  impose  the  boundary  conditions: 


p(0,t)  =  p(a,t)  =0 


The  validity  of  this  assumption  will  be  determined  by  the  success  of  our  fit.  The  solution 
expanded  in  a  Fourier  series  is  of  the  form'^: 


p(x,t)  =  Z  p„(t)  sin(n7tx/a)  n-^I,3,5,... 


where 


Pn(t)  =  [4ap</(nVD)]  sin(n7tx„/a)  (l-exp(-D(mi/a)h)) 


The  solution  of  the  diffusion  equation  after  the  discharge  is  turned  off  is  given  by  Eq. 


p(x,t-tj  =  I  p,(0  siii(imx/a)  exp(-D(nn/a)^(t-to)) 


n=l,3.5,... 


(8) 


where  is  the  time  the  discharge  is  turned  off,  and  p,(tj  is  given  by  Eq.  (7).  The  solid  line 
shown  in  Fig.  13  represents  a  semilogarithmic  fit  of  solution  (8).  Its  linearity  indicates  that  the 
first  term  of  the  expansion  of  the  solution  (8)  is  the  dominant  one  in  the  decay  of  the  atomic 
hydrogen  concentration  at  the  center  of  the  discharge.  The  dashed  line  which  resulted  from  a 
similar  fit  of  the  first  only  term  of  the  expansion  is  indicative  of  the  fact  that  higher  order 
terms  are  important  only  for  times  <  200  ps  while  they  contribute  less  than  5  percent  in  the 
decaying  concentration  at  times  >  200  ps  after  the  discharge  is  turned  off.  The  first  term  of 
the  series  relates  the  time  constant  r  of  the  exponentially  decaying  H-concentration  to  the 
diffusion  coefficient  D  :  T=(a/n)^/D.  We  thus  calculate  the  diffusion  coefficients  for  H-  and 
D-atoms  to  be  355±30  cmVs  and  220±15  cm  Vs  respectively.  Actually,  the  comparison  of  the 
diffusion  coefficients  of  these  atoms  provides  a  test  for  the  model  discussed  above,  since  the 
diffusion  coefficient  depends  inversely  on  the  square  root  of  the  mass  of  the  diffusing  particle. 
The  estimated  values  essentially  agree  with  the  expected  2^^  ratio  within  the  experimental 
uncertainty. 

The  steady-state  spatial  distribution  p„(x)  is  given  by  Eqs.  (6)  and  (7),  if  we  let  t->oo: 


•J 

♦ 


p„(x)  =  4apy(n^D)  2  (1/n^)  sin(nnXo/a)  sin(n7ix/a)  n=I,3,5,...  (9) 

We  fit  the  theoretical  spatial  distribution  to  the  experimental  data  with  parameters  D,  p„,  and 
x„.  The  best  fitting  curves  obtained  by  a  least  square  fit  program  are  shown  (along  with  the 
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data  used  for  the  fits)  in  Figs.  10,  11,  and  12.  We  estimate  the  position  of  the  sources  to  be 
at  3.4±0.2  mm  away  from  the  electrode  surfaces,  the  production  rate  to  be  (4.1±0.6)xl0'‘ 
atoms/cm^sec  and  7.8x10“  atoms/cm^sec  in  the  pulsed  and  the  CW  discharge  respectively 
(operated  at  5  W  and  10  W).  The  value  of  the  diffusion  coefficient  D  as  determined  in  these 
spatial  profile  fits  is  380±60  cmVs.  This  value  agrees  well  with  the  value  from  the  temporal 
decay  data  and  agrees  (within  the  uncertainty)  with  that  of  440  cm^/s  as  calculated  from 
Gilliland’s  semiempirical  equation."  This  equation  based  on  kinetic  theory  assumes  a  T^'^ 
dependence  for  the  diffusion  coefficient  and  it  is  evaluated  at  300  K  (we  assume  thermally 
equilibrated  H-atoms  at  room  temperature).  The  production  of  energetic  species  at  regions 
located  several  millimeters  away  from  the  surfaces,  has  been  indicated  previously  for  species 
such  as  F’  metastables, O  atoms,“  and  BCl  radicals"  at  high  frequencies.  It  was  shown 
in  the  last  work  that  the  production  of  radicals  is  affected  by  the  frequency  of  the  discharge. 
One  may  also  expect  that  the  pressure  affects  the  regions  of  production  of  energetic  species. 
Therefore  the  dependence  of  the  location  of  the  H-atom  sources  on  frequency  and  pressure 
should  be  the  subject  of  further  study. 

The  best  fitting  line  we  show  in  Fig.  10  represents  the  steady-state  spatial  solution  (9) 
for  totally  absorbing  surfaces  at  x=0,  x=a.  However,  the  experimental  concentration  gradient 
near  the  electrode  surfaces  is  smaller  than  that  calculated  from  Eq.  (9),  an  indication  that  the 
electrode  surfaces  are  partially  absorbing  with  a  loss  coefficient  y,  for  H-atoms: 

y  =  D(ap/5x)/[D(5p/^)+vp/4]  [  ,  (10) 
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where  D(5p/5x)  |  is  the  flux  of  H-atoms  that  are  absorbed  on  the  stainless  steel  electrode 
surface  at  x=0  and  vp/4  is  the  flux  of  H-atoms  of  mean  thermal  velocity  v  striking  on  the 
surface  without  sticking.  The  loss  coefficient  of  H-atoms  on  a  surface  can  be  determined  from 
semilogarithmic  fits  of  the  steady  state  spatial  profiles  near  the  electrode  surfaces,  as  those 
shown  in  Fig.  14.  The  fitting  curves  are  obtained  by  solving  the  diffusion  equation  for 
boundary  surfaces  with  a  loss  coefficient  y,=5  percent  for  hydrogen  and  deuterium  atoms.  The 
uncertainty  in  determining  y,  from  the  semilogarithmic  fits  is  ±1  percent.  The  loss  coefficients 
for  hydrogen  and  deuterium  are  equal,  as  one  would  expect  from  expression  (10)  for  y*  since 
the  mean  thermal  velocity  and  the  diffusion  coefficient  both  vary  as  m  '*.  The  result  further 
verifies  the  fact  that  the  balance  between  diffusion  and  atom  removal  at  the  surface  establishes 
the  observed  concentration  gradient  close  to  the  electrode  surfaces.  In  addition,  the  value 
D=355±20  cm  Vs  determined  by  the  temporal  decay  should  then  be  modified  to  account  for  the 
partially  absorbing  electrode  surface.  Numerical  solution  of  the  diffusion  equation  for  H-atoms 
confined  between  two  surfaces  with  a  5  percent  loss  coefficient  gives  a  coirected  value  of 
D=460±25  cmVs,  in  close  agreement  with  the  literature  value  D=440  cmVs. 

C.  H-ATOM  PROFILES  IN  LOADED  DISCHARGES 

i)  Experimental  observations 

In  Fig.  15,  we  show  hydrogen  atom  concentration  profiles  taken  1.5  mm  above  the 
grounded  bottom  electrode  while  the  reactor  cell  is  translated  horizontally  (radial  profiles).  The 
radial  profile  shown  in  Fig.  15(a)  is  taken  over  a  10x10  mm  GaAs  wafer.  We  observe  a 
significant  increase  in  hydrogen  atom  concentration  localized  over  the  wafer  surface.  A  radial 
profile  over  an  identical  stainless  steel  piece,  shown  in  Fig.  15(b),  proves  that  the  drastic 
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Figure  14.  Steady-state  concentrations  for  H-  and  D-atoms  at  small  distances  from  the 

stainless  steel  electrode  surface.  The  solid  lines  are  fits  obtained  for  surfaces  • 

4  with  a  loss  coefficient  y,«S%  for  hydrogen  and  deuterium  atoms.  Discharge 

conditions  the  same  as  in  Fig.  10. 
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Comparison  of  radial  profiles  of  H-atom  concentration  probed  1  mm  above  (a)  a  GaAs  wafer,  (b)  a  s( 
substrate,  and  (c)  a  gold  coated  GaAs  wafer.  A1  wafers  of  10  x  10x0.5  mm  size  are  centered  on  the  si 
ground  electrode  extending  from  -25  mm  to  +25  nun  of  the  x-axis. 


change  in  the  concentration  cannot  be  attributed  to  changes  in  the  bulk  properties  of  the 
discharge,  or  in  the  detection  geometry.  The  GaAs  surface  must  be  responsible  for  the  increase 
in  the  concentration,  as  we  further  verified  by  depositing  2-pm-thick  gold  on  the  GaAs.  A 
radial  profile  taken  over  the  gold  coated  GaAs  waf<^,  shown  in  Fig.  1 5(c),  indicates  a  reduction 
in  the  hydrogen  concentration  and  suggests  that  a  surface/atomic  hydrogen  interaction  must  be 
occurring. 

By  employing  a  pulsed  discharge,  we  studied  the  temporal  development  of  the  hydrogen 
concentrations  in  the  interelectrode  space.  In  Fig.  16  we  show  hydrogen  atom  concentrations 
in  a  discharge  loaded  with  a  15-cm^  GaAs  wafer  at  various  times  after  the  onset  of  the 
discharge.  Specifically,  the  spatial  concentrations  shown  were  measured  at  times  t=0. 1 5  ms, 
0.50  ms,  and  2.50  ms.  The  profile  at  t=0.I5  ms  indicates  that  the  production  of  hydrogen 
atoms  in  these  plasmas  occurs  in  localized  regions  in  close  proximity  to  the  electrodes.  In  time 
the  hydrogen  atoms  diffuse  from  their  sources  towards  the  center  and  the  electrode  surfaces. 
We  observe  that  the  spatial  profile  at  very  early  times  is  not  affected  by  the  presence  of  the 
GaAs  wafer;  however,  as  time  increases,  the  concentration  of  the  hydrogen  atoms  builds  up 
close  to  the  wafer  surface,  resulting  to  the  apparent  spatial  asymmetry  of  the  profile  taken  at 
time  t=2.5  ms. 

In  Fig.  17,  we  show  axial  profiles  of  hydrogen  atom  concentration  in  a  continuous 
discharge,  unloaded  and  loaded  with  Si  and  GaAs  wafers.  A  typical  axial  profile  in  a 
continuous  unloaded  hydrogen  discharge  is  shown  in  Fig.  17(a).  A  concentration  of  8x10'^ 
atoms/cm^  stays  relatively  constant  over  a  wide  region  in  the  interelectrode  space,  until  it 
exhibits  a  considerable  gradient  in  the  regions  0-4  mm  away  from  each  electrode  surface.  In 
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Figure  16.  Evolution  of  H-atom  spatial  distribution  in  a  pulsed  discharge  loaded  with  a  GaAs 
wafer  on  the  grounded  electrode;  experimental  data  points  and  simulated  results 
(solid  and  dashed  lines)  at  r=0. 15  ms,  r=O.SO  ms,  and  r-2.S0  ms.  The  location 
of  the  ground  and  powered  electrodes  is  a  0  and  19  mm,  respectively.  Discharge 
conditions:  pressure  of  3  Torr,  and  power  of  40  W. 
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H-atom  concentration  vs.  distance  from  tbe  grounded  electrode  in  a  continuous  discharge;  (a)  unloi 
(b)  a  Si  wafer,  and  (c)  a  GaAs  wafer  placed  on  tbe  grounded  electrode  (0-0.5  mm).  The  solid 
represent  fits  obtained  from  a  diffusion  model.  Discharge  conditions;  pressure  of  3  Torr,  and  po’ 


Fig.  17(b)  we  show  an  axial  profile  taken  under  identical  conditions  with  the  exception  that  a 
Si  wafer  is  placed  on  the  grounded  electrode,  covering  three-fourths  of  the  electrode  surface 
area.  We  observe  a  significant  qualitative  difference  from  the  previous  profile.  The  hydrogen 
concentration  appears  considerably  increased  in  the  region  from  the  center  of  the  discharge  to 
the  electrode  surface  that  bears  the  Si  wafer.  Similarly,  Fig.  17(c)  shows  a  spatially 
asymmetric  hydrogen  profile  in  a  discharge  loaded  widi  a  GaAs  wafer. 

ii)  Model 

We  explain  the  observed  behavior  of  the  H-concentration  with  an  extension  of  the 
model  for  an  unloaded  discharge,  according  to  which  these  profiles  are  determined  by  diffusion 
of  the  hydrogen  atoms  from  the  regions  where  they  are  produced  towards  the  surfaces  which 
confine  the  plasma.  The  model  assumes  two  point  sources  of  hydrogen  atoms  located  at  a 
distance  away  from  each  electrode  surface.  We  also  assume  that  u^t  differentiates  the 
behavior  of  the  atomic  concentration  at  the  presence  of  different  surfaces  in  the  discharge  is 
the  nature  of  the  surfaces  themselves.  Metallic  surfaces  are  known  to  enhance  the 
heterogeneous  reactivity  ^of  radicals  as  opposed  to  oxidized  semiconducting  surfaces.  Direct 
evidence  for  that  is  provided  in  Fig.  1 5,  where  H-atoms  are  consumed  near  metallic  surfaces 
(stainless  steel  and  gold)  rather  than  near  the  GaAs  surface.  Consequently,  we  assume  in  a 
first  approximation  that  the  stainless  steel  electrodes  represent  a  perfectly  absorbing  surface  for 
the  hydrogen  atoms  with  a  sticking  or  recombination  coefficient  close  to  1,  while  the 
semiconducting  surfaces  represent  totally  reflecting  surfaces.  We  again  solve  the  diffusion 
equation"  for  tfie  hydrogen  atom  density  p(x,t); 


ap/at  -  D  a^p/ax^  =  p„  (8(x-x„)+  8(x-a+x„)) 


(H) 

where  p^  is  the  production  rate  of  the  H-atoms,  D  is  the  diffusion  coefficient  for  H-atoms,  x„ 
is  the  distance  of  the  point  sources  from  each  electrode,  and  a  is  the  interelectrode  distance. 
For  the  unloaded  discharge  the  boundary  conditions  were  f>(0,t)-p(a,t)=0.  Above  a  wafer 
loaded  on  the  electrode  at  x=0,  we  set  the  boundary  conditions  ap/at(0,t)=0,  p(a,t)=0.  Then 
the  solution  of  Lq.  (II)  is: 

p(x,t)  =  2apy(7t^D)  Z  p„(x)  (l-exp(-(n+l/2)^(n/a)^Dt)]  n=0,I,2 .  (12) 

where 

p„(x)=(n+l/2)'^  [sin((n+l/2)rrXo/a)+(-l)"cos((n+l/2)7rxya)l  sin((n+l/2)7r(a-x)/a)  (13) 

We  discuss  the  validity  of  this  model  in  describing  the  data  in  Figs.  16  and  17.  The  solid 
lines  in  Fig  16  represent  the  solutions  (12)  of  the  diffusion  equation  at  t=0.15  ms  and  t=0.50 
ms.  For  both  fits  the  diffusion  coefficient  is  390  cmVs,  the  H-atom  sources  are  centered  at  3.4 
mm  away  from  each  electrode  surface,  and  the  prxxiuction  rate  of  the  hydrogen  atoms  is 
3.2x10”  atoms/cm^s.  The  values  of  the  parameters  D=390±60  cmVs  and  x„=3.4±0.2  mm  were 
also  used  to  successfully  fit  the  steady-state  concentration  profile  in  an  unloaded  discharge  (see 
Fig.  17(a)),  as  well  as  temporally  evolving  spatial  profiles  of  atomic  hydrogen  in  unloaded 
discharges.”  The  value  of  the  diffusion  coefficient  D  for  the  hydrogen  atoms  is  in  agreement 
(within  the  uncertainty)  with  the  value  of  440  cmVs  which  has  been  calculated  assuming  a  T’^ 


dependence  of  D."  The  fact  that  we  obtain  reasonably  good  fits  for  loaded  and  unloaded 
discharges  with  the  same  value  of  x,  leads  to  the  conclusion  that  tl^  presence  of  the 
semiconducting  surfaces  does  not  alter  the  location  of  die  sources  of  hydrogen  atoms,  at  least 
as  far  as  the  sensitivity  of  our  method  allows  us  to  measure. 

The  experimental  data  at  t=2.S0  ms  show  a  hydrogen  atom  concentration  substantially 
lower  than  what  the  model  predicts  for  a  perfectly  reflecting  surface.  A  partially  reflecting 
GaAs  surface  may  accommodate  this  discrepancy.  The  dashed  line  shown  in  Fig.  16  represents 
a  fit  which  uses  the  H-concentration  just  above  the  GaAs  surface  as  a  parameter.  We  found 
that  a  reduction  of  the  hydrogen  atom  concentration  adjacent  to  the  GaAs  surface  by  about  30 
percent  of  its  value  near  a  perfecUy  reflecting  surface  or  equivalendy  a  0.2  percent  sticking  or 
recombination  coefficient  for  the  GaAs  surface  provides  a  reasonably  good  fit  to  the  data.  This 
result  represents  a  recombination  coefficient  typical  for  native  oxides'*  inherent  on  the 
semiconducting  surfaces. 

The  solid  line  in  Fig.  17(a)  represents  the  solution  (3)  of  the  difliision  Eq.  (2)  at  steady 
state  ( t-Ko).  The  parameters  for  the  fit  are  0=390160  cmVs,  x„=3.4±0.2  mm  and  p,=  7.8x10'* 
atoms/cm^s.  The  production  rate  is  lower  than  that  in  the  pulsed  discharge  due  to  the 
smaller  power  we  used  for  the  continuous  discharge  (10  W)  as  opposed  to  a  power  of  40  W 
for  the  pulsed  discharge.  We  observe  that  the  model  is  in  good  agreement  with  the  data  for 
the  unloaded  discharge.  To  fit  the  data  taken  in  the  loaded  discharges  (Figs.  17(b),  17(c)),  we 
had  to  use  the  H-concentration  adjacent  to  the  Si  and  GaAs  surfaces  as  a  parameter,  as  we  did 
for  the  fit  in  Fig.  16(c).  The  dashed  lines  in  Figs.  17(b),  17(c)  represent  fits  with  the  same 
parameters  D,  and  a  decrease  of  the  concentration  of  H-atoms  adjacent  to  the  Si  and  GaAs 
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surfaces  by  about  40  percent.  The  agreement  between  the  simulations  and  the  experimental 
data  is  indicative  of  partially  reflecting  surfaces  or  equivalently  a  partial  loss  of  atomic  H  on 
the  surface.  The  consumption  of  H-atoms  may  further  rise  by  the  removal  of  native  oxides 
inherent  on  the  surfaces  in  discharges  at  lower  pressures  containing  heavier  etchant 
gases. Clearly,  further  investigation  is  required  to  illuminate  the  details  of  the  interaction 
between  the  hydrogen  atoms  and  the  semiconducting  surfaces. 

D.  LINEWIDTHS  AS  A  PROBE  OF  DISCHARGE  TEMPERATURE  AND  PURITY 

This  study  focussed  on  determining  the  temperature  of  H-atoms  in  the  RF  discharge. 
The  kinetic  energy  of  reactive  species  in  a  plasma  should  be  a  key  parameter  in  the  modeling 
of  reactor  phenomena  and  may  have  significant  effect  on  devices  or  materials  produced  by  the 
processing. 

Assuming  the  H-atoms  are  produced  by  the  discharge  itself,  conventional  wisdom 
suggests  that  such  a  relatively  long-lived  species  will  be  equilibrated  with  the  overall  heavy 
particle  temperature.  However,  previous  work^'  utilizing  plasma  induced  emission  has 
determined  that  a  large  fiction  of  the  H-atoms  in  excited  electronic  states  are  very  hot  indeed, 
with  temperatures  in  excess  of  10*  K.  The  present  experiments  provide  the  first  comparable 
measurement  for  ground  state  H-atoms,  obviously  the  numerically  predominant  form  present 
in  the  plasma. 

In  the  emission  work,  the  mechanism  heating  the  excited  state  atoms  is  charge  exchange 
with  fast  moving  (hot)  hydrogen  ions,  H*,  Hj*,  Hj*,  accelerated  by  the  electric  fields  in  the 
discharge.  A  similar  ion  exchange  mechanism  almost  certainly  exists  for  ground  state  H-atoms. 
The  question  is  whether  their  somewhat  longer  lifetime  is  sufficient  to  equilibrate  the  vast 
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majority  of  them  to  ambient  temperature. 

Another  mechanism  that  will  produce  hot  H-atoms  is  Hj  electron  impact  dissociation 
process  in  the  plasma  which  is  well  known^'”'^*'^’  to  produce  atoms  with  translational 
energies  of  order  10^  K.  Again,  the  question  is  whether  these  hot  atoms  are  equilibrated  to 
ambient  temperature. 

A  final  process  which  can  produce  hot  H-atoms  is  the  TALIF  diagnostic  itself 
Photodissociation  of  a  molecule,  like  electron  impact  dissociation,  can  produce  very  hot  atoms. 
Since  photodissociation  cross  sections  of  ground  state  by  20S-nm  photons  are  small, 
photodissociation  by  the  diagnostic  laser  is  a  negligible  process  in  a  pure  plasma.  However, 
many  hydrogen  containing  molecules,  e.g.  NH3,  hydrocarbons  from  pump  oil  and  elsewhere, 
etc.,  have  bond  energies  less  than  that  of  the  20S-nm  diagnostic  photons  and  appreciable  cross 
sections.  Should  those  molecules  be  present,  either  as  feedstock  gases,  contaminants,  or  via 
plasma  synthesis,  then  photodissociation  leading  to  hot  H-atoms  is  quite  likely.  Moreover, 
these  H  atoms  will  likely  be  detected  before  they  have  any  opportunity  to  undergo  equilibrating 
collisions,  since  they  may  be  observed  by  subsequent  205-nm  photons  in  the  stune  laser  pulse 
whose  duration  is  ~5  nsec,  shorter  than  the  lifetimes  of  the  hot  excited  state  H-atoms 
previously  observed  to  be  produced  by  ion  exchange.  Thus  H-atoms  produced  photolytically 
will  bear  a  spectral  profile  signature  characteristic  of  their  precursor  molecule. 

During  our  work  we  have  analyzed  H-atom  TALIF  signals  produced  both  in  an  RF  H^ 
plasma  and  downstream  from  an  Hj  microwave  plasma.  We  have  also  examined  RF  plasmas 
containing  other  processing  gases  such  as  ammonia,  acetylene  and  hydrogen  sulfide,  and  these 


results  will  be  discussed. 


i)  Theory 

Our  procedure  for  determining  H-atom  ten:^)eratures  is  to  measure  the  TALIF  linewidth. 
Under  ideal  conditions,  the  experimental  line  profile  will  be  entirely  determined  by  Doppler 
broadening.  In  this  case,  the  full  width  at  half  height,  Avq,  is  related  to  the  atomic  velocity  and 


temperature  by^‘ 


AVo=2vo(v/c)=(2v</c)|^^^j''*=7.16xlO-’(77Af)‘^v 


where 


V  =  average  velocity  of  the  atom  or  molecule 
c  =  speed  of  light 
v^  =  transition  frequency 
k  =  Boltzmann  constant 
m  =  mass  of  the  atom  or  molecule 
T  =  temperature  of  the  atom  or  molecule 
M  =  molecular  weight  in  amu  of  the  atom  or  molecule 
To  give  one  an  appreciation  of  the  magnitudes  involved,  it  is  useful  to  evaluate  Eq.(14)  for  H- 
atoms  at  ambient  temperature  (300K).  The  experimentally  assessed  TALIF  3  ^S(^D)  o  1 
H-atom  transition  occurs  at  v„  =  97492  cm  '.  Thus  Av,  =  1.21  cm  ’  for  ambient  H-atoms. 
Considering  the  fact  that  we  combine  6  dye  laser  photons  to  drive  this  transition  implies  that 
the  maximrim  width  of  our  dye  laser  ^ouid  be  ^  0.1  cm  '  to  produce  a  predominantly  Doppler 

broadened  line.  Additionally,  we  note  that  because  the  molecular  weight  occurs  in  the 


denominator  of  Eq.(I4),  the  requirement  on  the  dye  laser  linewidth  will  be  much  more  stringent 
to  obtain  temperatures  of  other  reactive  atoms. 

it)  H2  Discharge 

To  establish  a  reference  for  measuring  ambient  TALIF  linewidths,  microwave  generated 
hydrogen  atoms  were  detected  after  flowing  down  a  few  centimeters  of  tubing.  The  microwave 
discharge,  well  removed  from  the  region  of  fluorescence  capture,  does  not  contribute  the  bright 
PIE  present  in  the  RF  plasma  cell;  therefore,  a  continuous  wave  photomultiplier  tube  was  used 
to  initially  record  the  H-atom  resonance.  This  was  a  test  to  ensure  that  the  response  time  of 
the  gated  tube  did  not  bias  our  linewidth  measurements.  Using  this  configuration,  a  linewidth 
of  1.34  ±  0.07  cm  '  (fwhm)  was  measured  for  low  laser  powers  (<  30  pJ).  Treating  the 
lineshape  as  Gaussian  and  deconvoluting  a  laser  linewidth  of  0.42  cm  '  (fwhm,  3  photons  to 
produce  205-nm  x  2-photon  process;  6  x  0.07  cm'')  results  in  an  inherent  Doppler  width  of 
1.27  cm  ',  which  corresponds  to  a  temperature  from  Eq.(14)  of  330  ±  35K..  Although  this  is 
slightly  hotter  than  room  temperature,  it  agrees  within  experimental  uncertainty.  A  small 
contribution  from  laser  saturation  broadening  may  also  add  slightly  to  linewidth. 

Linewidths  measured  with  the  GPMT  under  the  same  conditions  were  1.47  ±  0.08  cm  '. 
This  value  was  very  sensitive  to  integration  parameters,  for  example,  gate  width  and  gate 
position.  A  gatewidth  of  75  ns  and  gate  position  centered  on  the  resonance  produced  the 
narrowest  measurements.  The  slight  discrepancy  between  the  CWPMT  and  the  GPMT  is 
attributed  to  the  response  of  the  gated  tube.  The  GPMT  has  a  tendency  to  easily  saturate, 
giving  the  ^pearance  of  a  slightly  broadened  line.  The  value  of  1.47  cm''  was  taken  to 
correspond  to  the  ambient  temperature  inclusive  of  instrumental  broadening  for  the  GPMT. 


The  simple  Gaussian  deconvolution  leads  to  an  effective  instnimental  contribution  of  0.6  cm 

Determination  of  the  temperature  for  the  RF  discharge  H-atoms  followed  the  same 
procedure  used  to  record  the  narrowest  H-atom  signals  produced  from  the  microwave 
discharge.  For  low  discharge  powers  and  low  discharge  pressures  (»0.5  torr,  <12  W)  a 
minimum  linewidth  of  1.64  ±  0.2  cm  '  is  obtained,  corresponding  to  a  temperature  of  440  ± 
HOK  after  instrumental  and  laser  linewidth  contributions  are  considered.  A  representative  H- 
atom  spectrum  is  shown  in  Fig.  18. 

The  relatively  large  uncertainty  in  the  temperature  measurement  in  the  RF  discharge, 
caused  by  PIE  and  by  using  low  laser  powers  to  minimize  saturation,  can  arguably  equate  the 
microwave  and  RF  temperature;  however,  other  evidence  suggests  that  the  RF  environment  is 
hotter  than  ambient  under  most  conditions.  When  the  rf  power  or  the  pressure  is  increased 
from  the  minimum,  the  linewidth  and  hence  the  temperature  also  increased  to  a  maximum  of 
over  lOOOK  (see  Fig.  19).  Deuterium  data  were  also  recorded  and  showed  a  linewidth  2  *  times 
nanower  than  H-atoms,  ensuring  that  the  increase  in  width  is  related  to  Doppler  broadening. 

Data  for  microwave  and  RF  discharge  produced  H-atoms  were  also  recorded  employing 
the  dye  laser  without  the  intracavity  etalon.  These  grating  scans  served  primarily  as  survey 
data,  but  reflect  the  same  trends  for  the  pressure  and  RF  power  dependence  on  linewidth,  when 
the  much  larger  laser  linewidth  contribution  is  deconvoluted. 

Such  high  temperatures  for  the  H-atoms  are  somewhat  surprising,  especially  considering 
that  previous  measurements  for  temperatures  of  BCl”  and  O*  atoms  in  Rf  plasmas  are  300K. 
One  would  expect  that  the  neutral  species  should  be  equilibrated  to  the  ambient  temperature 


of  the  reactor  walls. 


igure  18.  Etalon  narrowed  TALIF  scan  of  the  H  atom  3*S(^D)«*1^(S(97,492.3  cm'')  transition  from  atoms  generated  in  an  RF 
discharge.  Actual  data  pcints  and  a  Gaussian  fit  to  the  data  are  presented. 


LInewidth  (cirr^) 


Pressure  (torr) 


Figure  19.  Comparison  of  linewkith  dependence  on  H2  pressure  (recorded  at  12  W  RF 

power)  and  on  discharge  power  (recorded  at  an  Hj  pressure  of  0.42  Torr).  Also  i 

shown  is  the  temperature  scale  corresponding  to  the  observed  linewidths. 
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The  only  possible  broadening  mechanisms  inherent  to  the  RF  discharge  and  exclusive 
of  the  microwave  discharge  would  appear  to  be  Stark  broadening  and  E>oppler  broadening. 
Both  mechanisms  also  appear  consistent  with  the  observed  pressure  and  discharge  power 
dependence.  However,  Stark  broadening  appears  unlikely  because  of  the  fact  that  D  linewidths 
are  narrower  by  2.5.  Additionally,  Stark  broadening  calculations  based  on  the  quasistatic 
microfield  approximation"  dictate  that  the  electron  density  in  the  plasma  must  be 
approximately  3  orders  of  magnitude  higher  than  the  10"  electrons/cm  ^  maintained  under 
typical  operating  conditions  to  produce  the  observed  broadening. 

The  only  plausible  explanation  is  Doppler  broadening.  We  propose  that  H-atoms  are 
being  warmed  by  the  charge  exchange  with  hot  ionic  species  in  the  glow.  Donnelly  et  al.^^ 
have  shown  that  by  increasing  the  power  to  the  discharge,  the  number  of  ions  increases  and 
then  levels  out  or  saturates  at  higher  power.  Figure  19  clearly  shows  a  similar  phenomena  in 
the  linewidth  vs.  power  dependence.  The  overall  temperature  varies  from  approximately  440 
to  900K  from  lower  to  higher  discharge  powers.  Increasing  the  pressure  also  increases  the  ion 
density  up  to  a  point  and  then  levels  off  producing  similar  results  for  the  linewidth  as  also 
shown  in  Fig.  19. 

Hi)  Ammonia 

In  contrast  to  the  above  results,  translationally  hotter  hydrogen  atoms  having  a 
linewidth  (see  Fig.  20)  of  approximately  5  cm  '  (fwhm)  corresponding  to  a  temperature  of 
*5000  K,  were  initially  observed  while  running  a  pure  Hj  RF  discharge.  However,  this 
discharge  was  later  found  to  be  contaminated  by  a  small  leak  in  the  vacuum  system. 
Remarkably,  the  TALIF  signal  remained  for  long  periods  (1  or  2  h)  after  the  discharge  was 
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Comparison  of  the  linewidth  of  RF  generated  H-atoms  (inner  trace)  and  atoms  generated  from  the  photodissociation 
of  ammonia.  Actual  data  points  and  Gaussian  fits  are  shown. 


extinguished.  It  is  postulated  that  H-atoms  react  with  Nj  or  N  atoms  to  form  NH3,  which  is 
subsequently  photolyzed  by  one  20S-nm  photon  with  the  resulting  H-atom  detected  by  TALIF. 
The  H-atom  would  carry  translational  energy  produced  by  the  excess  energy  of  the  photolysis 
laser  relative  to  the  N-H  bond  dissociation  energy,  Dg.  Photolysis  and  excitation  occur  during 
the  same  pulse  (S  ns),  eliminating  the  possibility  of  thermalizing  collisions. 

Earlier  work  by  Back  and  Koda^*  shows  that  dissociation  of  NH3  is  easily 
accomplished  by  absorption  of  a  photon  of  20S'nm  light  to  a  predissociative  state.  Results  of 
Uyama  and  Matsumoto”  show  that  NH3  is  synthesized  in  microwave  discharges  of  Hj/Nj 
mixtures  by  formation  of  NH  radicals  that  further  react  with  either  H  or  Hj. 

Confirmation  of  the  theory  of  NH3  synthesis  was  established  by  comparison  of  TALIF 
linewidths  for  Hj/N,  mixtures  with  TALIF  linewidths  for  NH3  held  static  in  the  cell  with  the 
discharge  off.  Results  for  TALIF  linewidths  from  hydrogen-nitrogen  RF  mixtures  and  pure 
ammonia  (given  in  Table  1)  are  identical  within  experimental  error.  Additional  evidence  is 
shown  in  Fig.  21,  which  gives  a  log  plot  of  H-atom  TALIF  intensity  versus  laser  power.  A 
slope  of  two  is  predicted  for  TALIF  of  discharge  generated  H-atoms,  while  a  slope  of  three  is 
expected  for  the  three-photon  process  of  photodissociation  followed  by  two-photon  absorption. 
As  Fig.  21  shows,  the  results  from  the  NH3  experiments  indicate  a  slope  close  to  3,  except  at 
highest  laser  powers  where  saturation  occurs.  It  is  clear  from  these  data  that  translationally  hot 
H-atoms  are  produced  by  photodissociation  of  NH3. 

The  translational  temperature  of  hot  atoms  produced  photolytically  provides  a  means 
of  identifying  the  precursor  species,  via  its  bond  dissociation  energy.  Segall  et  al.^ 
demonstrated  a  technique  for  measurement  of  bond  dissociation  energies  from  lineshapes  of 
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TABLE  1  Processing  Gas  Bond  Energies  and  Experimental  Linewidths 


LinewWtlifca') 


5.01(10) 


5.96(14) 


4.48(8) 


“values  taken  from  ref.  29. 
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H-atoms  photolyzed  from  acetylene.  This  technique  yields  a  controversial  result,^'  but  is  at 
worst  a  good  iq)proximation  (±  10%).  The  observed  lineshape  is  followed  out  into  the  wings 
to  determine  maximum  kinetic  energy  of  the  H-atom.  The  maximum  velocity,  vJJ  of  the  H- 


atoms  contributing  to  the  Doppler  profile  can  be  easily  related^^  to  the  bond  energy.  Dp, 

Dfl  ~  E205  -  Ejr  =  E205  -  SEh  =  Ejoj  -  1/2  Sm||(v^)  ,  (15) 

where  Ep  is  the  energy  of  ail  fragments,  E^oj  is  the  photolysis  energy,  and  E^  is  the  maximum 
energy  of  the  hydrogen  atom  fragment.  The  numerical  factor  S  assures  conservation  of 
momentum  in  the  photodissociation  process.  It  is  entirely  determined  by  the  masses  of  the 
fragments;  for  NH,,  CjHj  and  HjS,  it  is  respectively  1.06,  1.08  and  1.03.  We  have  adopted 
the  procedure  of  fitting  the  H  TALIF  lineshape  to  a  Gaussian  profile  and  determining  E„  at  the 
2.5  sigma  point  on  the  curve.  We  were  able  to  derive  a  bond  dissociation  energy  for  H-NHj 
of  109  ±  11  kcal/mole  which  agrees  well  with  the  accepted  value  as  shown  in  Table  1. 

ivj  Other  Processing  Gas  Plasmas 

The  quality  of  the  determined  bond  dissociation  energy  for  ammonia  led  to  studies  of 
other  processing  gases  to  establish  whether  they  may  be  identified  by  lineshape  analysis. 
TALIF  signals  for  hydrogen  liberated  from  HjS  and  acetylene  were  easily  obtained  and 
linewidths  and  bond  dissociation  energies  are  also  presented  in  Table  1.  Fig.  22  shows  the 
contrast  between  the  lineshtqies  of  HjCj  and  HjS  generated  H-atoms.  From  Table  1  and  Fig. 
22,  it  is  clear  that  the  TALIF  diagnostic  can  identify  the  precursor  molecule. 

Attempts  were  also  made  to  observe  TALIF  signals  from  photodissociation  of  methane 
but  were  unsuccessful.  This  is  presumably  because  methane  does  not  possess  sufficient 
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oscillator  strength  for  absorption  of  a  205-nm  photon,  despite  the  photon  having  sufficient 
energy  to  rupture  the  C-H  bond.””  The  same  results  were  also  observed  with  Hj, 
presumably  for  the  same  reason.  This  conclusion  is  strengthened  by  the  almost  perfect 
correlation  shown  in  Fig.  21  between  H-atom  signal  and  squared  laser  power  for  RF  generated 
H  atoms  in  an  H2  discharge. 

While  not  a  perfect  method  for  identifying  species  in  plasma  environments,  the 
photodissociation  method  should  prove  helpful  in  this  respect.  We  found  that  dirty  tubing, 
dirty  valves,  and  even  the  slightest  leaks  in  our  vacuum  system,  resulted  in  broadening  of  our 
TALIF  signals.  Thus  in  a  processing  reactor  with  TALIF  diagnostic,  the  width  of  the  H-atom 
signal  can  provide  a  very  sensitive  indicator  of  its  proper  functioning.  We  have  certainly  found 
it  to  be  so  with  our  reactor. 
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SECTION  4 


CONCLUSIONS 

We  have  successfully  developed  FALIF  as  a  highly  spatially  and  temporally  resolved 
diagnostic  of  H-atom  relative  and  absolute  concentrations  in  cold,  reactive  plasmas,  such  as  are 
commonly  used  for  processing  of  semiconductors  and  other  materials  {4)plications.  We  have 
developed  a  model  which  accurately  accounts  for  variations  of  H-atom  concentrations  both  in 
sp>ace  and  time  for  plasmas  with  and  without  semiconducting  wafers  present.  We  have  further 
shown  how  H-atom  TALIF  can  be  used  as  a  sensitive  probe  of  the  temperature  and  purity  of 
the  plasma. 

These  experiments  have  demonstrated  the  usefulness  of  the  TALIF  diagnostic  both  as 
a  means  for  understanding  the  fundamentals  of  the  discharge  and  for  monitoring  plasma 
processes.  It  is  clear  that  the  TALIF  techniques  can  be  extended  to  monitor  other  atomic 
species,  O,  N,  Cl,  etc.,  important  in  plasma  processing.  Furthermore,  the  high  spatial  resolution 
should  be  very  useful  for  probing  near-surface  interactions.  To  further  generalize  the 
usefulness  of  the  technique,  it  will  be  important  to  extend  measurements  to  higher  pressure 
plasmas.  To  do  this,  reliable  measurements  of  quenching  crosssections  for  the  fluorescence 
must  be  made. 
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